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Abstract 

A dynamical system with discrete time is studied by means of algebraic 
geometry. The system admits a reduction that is interpreted as a classical 
field theory in 2+1-dimensional wholly discrete space-time. The integrals 
of motion of a particular case of the reduced system are shown to coincide, 
in essence, with the statistical sum of the well-known (inhomogeneous) 2- 
dimensional dimer model (the statistical sum is here a function of two 
parameters). Possible generalizations of the system are examined. 

Vacuum curves and vacuum vectors are algebro-geometrical objects that 
have arisen in the theory of the quantum Yang — Baxter equation. They seem 
to have their origin in R. Baxter's works [Q, ||, and their general definition was 
formulated by I. M. Krichever Q. Baxter, and also Takhtajan and Faddeev 
[Q used the vacuum vectors to obtain a generalization of the Bethe ansatz for 
the XYZ spin model — an integrable one-dimensional quantum field theory (and 
also for the eight-vertex model of two-dimensional statistical physics) . Krichever 
has applied the vacuum curves to classification of the solutions of the quantum 
Yang — Baxter equation in the tensor product of 2-dimensional vector spaces. 
Then, the author of this paper has found some further applications of the vac- 
uum curves and vacuum vectors. In the paper ^] (see also H), new solutions of 
the quantum Yang — Baxter equation were constructed for the first time. They 
correspond to what is now known as Chiral Potts model. In jl0| O [)), the 
degeneracies of the spectrum of the XXZ quantum chain hamiltonian were ex- 
amined by means of the vacuum curves, and in []l2| , ^ the studying of vacuum 
vector bundles has resulted in the construction of solutions to the tetrahedron 
equation with commuting spin variables on the links. 

Here, I try to demonstrate that the vacuum curves may be useful for studying 
the classical (not quantum) field theory models as well. A difference equation 
on the 2+1-dimensional cubic lattice is presented, for which the solution to the 
Cauchy problem is constructed, at least in principle, through a rather simple 
scheme. The evolution is of hyperbolic nature, i.e. the "perturbations" propa- 
gate not faster than fixed speed. The interesting feature is that, in a particular 
"scalar" case, the model reveals a quite natural connection with the well-known 
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dimer model of statistical physics. The statistical sum of this latter model, 
which depends here on two parameters (and the model itself is, of course, inho- 
mogeneous), is the integral of motion for any values of these parameters. 

This field theory comes as a "reduction" of some very simply described 
"non-local" dynamical system. On the other hand, generalizations of this latter 
system are constructed in this paper, and I use a discrete analog of Lax pair for 
this purpose. 

Acknoledgements. I owe to A. B. Shabat the idea of "local reduction" 
(see Section ||). L. D. Faddeev informed me of the paper fl3|| . I would like to 
express my gratitude to them. 



1 Definition of the dynamical system. Gauge 
invariance 

Let 

L-f A B 
\ C D 

be a block matrix, A, . . . D being nx n matrices consisting of complex numbers. 
Consider the following two operations: construction of the inverse matrix 



and the block transposing 



CD) \ B D 

Now let a (birational) mapping / be a composition of these two operations: 

f(L) = (L- 1 )'. (1) 

Let us introduce the discrete integer-valued time r, and let the matrix L depend 
on r so that 

L(r + l)=/(i(r)). (2) 

This "dynamical system" has been already mentioned in literature p3[ . In 
the present paper, the integrability of this system is demonstrated, assuming 
that the "motion" is considered up to a "gauge transformation" (see below). 

Let G and H be non-degenerate n x n matrices. The gauge transformation 
of the matrix L is the following transformation of its blocks: 

A -> GAH, ... D ^ GDH. (3) 

Two matrices L and L' connected by the transformation (^) will be called gauge 
equivalent. It is clear that if L(t) and L'(t) belong to the same class of gauge 
equivalence, the matrices L(t + 1) and L'(t + 1) also do so. Thus, dynamics (^|) 
induces a dynamics on the set of classes of gauge invariance. 
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2 Vacuum curves and vacuum vectors 



It turns out that the dinamics @ preserves the so-called vacuum curve T of 
the operator L (the bases being fixed, we make no difference between a linear 
operator and its matrix). To be exact, T remains unchanged under the trans- 
formation / o /, and undergoes a simple transformation under /. The curve T 
together with the class of linear equivalence of the pole divisor of the vacuum 
vectors (see below) determines the matrix L up to a gauge transformation. The 
set of those classes of linear equivalence is isomorphic to a complex torus — the 
Jacobian of the curve T. The dynamics (^) linearizes on the Jacobian, i.e. the 
transformation / corresponds to a constant shift on the torus. Now, let us 
discuss these facts in detail. 

The vacuum curve of the operator L is an algebraic curve in the space C 2 of 
two variables u,v. Here are two equivalent definitions of it ||. 

Definition 1 Consider the relation 

L(U®X) = V®Y, (4) 



wherein 



are two-dimensional vectors, X and Y are n-dimensional vectors. For a 
generic matrix L, the non-zero solutions (U, V, X, Y) of the relation (Qj are 
parametrized, up to a scalar factor in X and Y , by points of an algebraic curve 
r of genus g = (n — l) 2 given by an equation of the form 

P(u,v) = 0, (5) 
P(u,v) being a polynomial of degree n in each variable, i.e. 

n 

P{u,v) = ^v k . (6) 

];k=l 

r is called the vacuum curve of the operator L. 

Definition 2 The vacuum curve of the operator L is the curve T in C 2 given 
by the equation 

P{u, v) = det(V ± LU) = det(uA + B- uvC - vD) = 0, (7) 

where 

V ± = (l,-v). 
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Let us denote the points of the vacuum curve by the letter z = (u, v) G 
T. Then U = U(z) and V — V(z) are meromorphic vectors on T with the 
pole divisors D\j and Dy of degree n, while X = X(z) and Y = Y(z), if 
normalized by, e.g., the condition that their nth coordinates equal unity, become 
meromorphic vectors with pole divisors Dx and Dy of degree n 2 — n p| . Under 
this normalization, a meromorphic scalar factor h(z) must be added into (Q): 

L{U{z)®X{z)) = h{z)V{z)®Y{z). (8) 

The linear equivalence of divisors 

D v + D x ~ Dy + Dy 

holds and is provided by the function h{z) in the sense that h(z) has its poles 
in the points of Djj + Dx and zeros in the points of Dy + Dy . 

As is shown in the paper ||, the vacuum curve equation P(u, v) = and the 
class of linear equivalence of divisor Dx or Dy determine a generic matrix L 
to within a gauge transformation, and vice versa, the gauge transformations do 
not change the vacuum curve and the classes of linear equivalence of divisors. 
In other words, the correspondence 

(class of gauge equivalence of L) <-> (r, the class of Dx) 

is a birational isomorphism. 

We will call X[z) the vacuum vector and Y[z) the covacuum vector in the 
point z of the curve T. X{z) = X(u,v) generates the (one-dimensional) kernel 
of the matrix 

uA + B - uvC - vD. (9) 

The Definition [l] allows one to trace what happens with the vacuum curve 
and vacuum vactors under the transformation L — > L^ 1 , while the Definition || 
allows one to trace what happens under the transformation L — * L . Namely, 
it is seen from the relation 

L-^Viz) <g> Y(z)) = hiz^Uiz) ® X{z) 

that the vacuum curve equation for the matrix L _1 is 

P(v,u) = 0, 

while its vacuum vector in the point (v, u) coincides with the covacuum vector 
of the initial matrix L. As for the block transposing, the vacuum curve equation 
for the matrix L 1 

det{uA + C - uvB - vD) = 

may be rewritten as 

u n v n &et{v~ 1 A -B + u- l v~ l C - u~ l D) = 0, 
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i.e. 

u n v n P(-v~\-u- 1 ) = 0. 

The vacuum vector of the matrix L l in the point (— v , — u^ 1 ) of its vacuum 
curve coincides with the vacuum vector X(u, v) of the matrix L. 

Combining these considerations, one finds out that the vacuum curve T of 
the matrix (X" 1 )* is given by equation 

u n v n P(-u- 1 ,-v- 1 ) = 0, 

while the vacuum vector X(— u v^ 1 ) coincides with the vector Y(u,v) of 
the matrix L. 

Identifying the curves T and T by means of the isomorphism 

(u,v) <-> (-u~\ -v' 1 ), 

one sees that 

D x ~ Dy ~ Dx + Du — Dv, 

which means that, in essence, the transformation ([!]) results in adding a fixed 
element of the Picard group, namely the equivalence class of the divisor Djj — 
Dy, to the pole divisor Dx of the vacuum vectors. It is clear also that after 
two transformations one returns to the initial curve: 

f = r. 

3 Reduction to evolution equation in the 2+1- 
dimensional space-time 

The dynamical system of the previous section admits an interesting reduction, 
i.e. some special choice of the matrices A, . . . D that is in agreement with the 
evolution. In this section, it will be convenient to treat the matrices A, ... D as 
linear operators acting from the linear space Hi into the linear space H2 (of the 
same finite dimension). This being the situation at the moment r, the operators 
act, of course, from Ti.2 into TL\ at the moment r + 1, and so on. 

Let each of the spaces Hi, TL2 be a direct sum of lm/2 identical subspaces 
of dimension d, where I, m are even numbers. Let us imagine these subspaces as 
situated at the vertices of the square lattice on the torus of the sizes I x m (which 
will mean the periodic boundary conditions in both discrete space variables). 
Let the subspaces be arranged in checkerboard fashion, as in Fig. ll where the 
empty circles correspond to subspaces of the space Hi, while the filled circles 
correspond to those of the space Hi- 

Let then the operators A, ... D be such that the image of each of the men- 
tioned (i-dimensional subspaces with respect to, say, operator A lies in the d- 
dimensional subspace of H2 at which points the arrow marked "A" that links 



5 




G 



However, the formulae (|10| 
After the transformation 



these two subspaces (Fig. |IJ). Analogously, the restrictions on B, C, D are de- 
picted in Fig. |l] (see also formula ( |20| ) for non-degenerate A, . . . D). Thus, to 
each link of the lattice a d x d matrix is attached that is a block of one of the 
"large" matrices A, . . . D. Let us shade half of the squares of the lattice in a 
checkerboard way, as in Fig. [I]. One can verify that the evolution of the system 
may be described as follows. 

At the first step, each of the four dxd matrices that correspond to the arrows 
surrounding each shaded square is transformed into a matrix expressed through 
just these four matrices. This goes according to the following formulae, in which 
the dxd blocks are somewhat freely denoted by the same letters A, . . . D as the 
"large" matrices: 

A — ► (A-BD^C)- 1 , (10) 

B — > (B- AC^D)- 1 , (11) 

C — ► (C-DB^A)- 1 , (12) 

D — > (D - CA^B)- 1 . (13) 

apply equally to the "large" matrices. 
131), all the arrows reverse, and at the second 
step the non-shaded squares are engaged in the same way according to the 
same formulae (|l0|-|l3|) . Then everything is repeated. Thus, the evolution is of 
hyperbolic nature: each local perturbation spreads not faster than one unit of 
length per unit of time. 

Let us clarify the symmetries of vacuum curves and divisors Dx in this 
"reduced" model. Let us introduce two integer-valued coordinates £, 77 for the 
vertices of the lattice, so that £ increases by 1 in passing from a vertex one 
step to the right, and 77 increases by 1 in passing one step upwards. £ and rj are 
defined modulo / and m respectively. A d-dimcntional subspace of "Hi or H.2 will 
be denoted Ti.^ v if it corresponds to a vertex with coordinates £, r\. Consider a 
linear transformation in spaces TL\ and TL2 consisting in multiplying the vectors 
of each subspace Tt^ by uf, oji being a fixed primitive root of the l-th degree 
of unity: 

uj[ = 1. 

This corresponds to the following transformation of the operators A, . . . D (from 
now on we speak of each of these operators "as a whole", not of their blocks): 

A -> uiA, B -> B, C -> C, D -> uj^D. (14) 

Consider also another linear transformation in Ti\ and Ti.2, consisting in multi- 
plying the vectors of each subspace Ti^ v by io\ , u>2 being a fixed primitive root 
of the m-th degree of unity: 

^ = 1. 

This corresponds to the following transformation: 

A -> A, B -> lj 2 B, C -> uj^C, D -> D. (15) 
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The vacuum curve of the operator L, which is given by equation (Q) 
P(u, v) = det(uA + B - uvC - vD) = 0, 



must be invariant under the transformations (J14|), (15). This leads to the invari- 
ance of the polynomial P(u, v) with respect to the following transformations g\ 
and g2- 

9i(u,v) = (uj 1 u,uJi 1 v), (16) 

92(u,v) = (uj 2 ; 1 u,uj2 1 v). (17) 

This invariance, then, leads to the following statement: only those coefficients 
a,jk are non-zero in the vacuum curve equation (see (g), @) for the "reduced" 
model, for which 

j-k = 0(mod /), 1 
j + k = 0(mod m). J 



(18) 



As for the divisor Dx, let us recall that it consists of such points in the curve 
T in which vanishes the last coordinate of the vector X (see [|| ) , the latter being 
an eigenvector of the matrix (^) with zero eigenvalue: 

(uA + B-uvC-vD)X(u,v) = 0. (19) 

This immediately leads to the conclusion: the divisor Dx is invariant with 
respect to the transformations (|l6|, |l7| ). 

Under some additional condition, the inverse statement also holds: if the 
curve r and divisor Dx are invariant under the transformations ([l6|), ([l7|), then 
the corresponding L-operator comes from a "reduced" model described in this 
section. For instance, this is true if 1/2 and m/2 are relatively prime numbers. 
If these numbers are not relatively prime, some conditions are to be imposed 
on the divisor Dx- To avoid going into details of this latter case, let us not 
consider it here. 

' A B 



Thus, let an operator L — I ^ ^ J be given, A, . . . D being n x n matrices, 

n = (lm/2)d, I and m even, and 1/2 and m/2 being relatively prime. Let the 
vacuum curve Y of the operator L and the divisor Dx be invariant under the 
action of the group Q generated by its elements gi, 32 ( |l6| , p"7j ), u>i and 0J2 being 
primitive roots of degrees I and m of unity. Then the linear space in which 
operators A, . . . D act decomposes into a direct sum of lm/2 d-dimensional 
subspaces Ti^, £ and rj being integers modulo I and m respectively and such 
that £ + 77 is an even number, and the following equalities between the images 
of these subspaces hold (in a "generic" case of non-degenerate A, . . . D): 

AH^-i,n+i = BTL^ri = CH^ :V +2 = DH^+i^+i- (20) 
The equalities (|20| ) mean exactly that one is in the situation of Fig. [|. 



8 



Let us prove the above statements. First, the natural projection from the 
curve r to its factor Y/Q has no branch points (here the fact that 1/2 and m/2 
are relatively prime is used to demonstrate that ramification does not occur 
when uoru equals zero or infinity). Thus, the n-dimensional linear space of 
meromorphic functions x(z) = x(u, v) whose pole divisor is Dx decomposes into 
a direct sum of subspaces of equal dimensions corresponding to the characters 
of (commutative) group Q. Each of these subspaces consists of functions x(z) 
satisfying relations 

x{gz) = x^{g)x{z), 

the character x^ being a scalar factor 

Xivi.9) = u\ a ujf, 

where 

9 = 9x9l 

The equality g^g™ = 1 means that £ + rj must be an even number. 

The components of the vector X(z) are exactly the functions x(z). In an 
appropriate basis, d components correspond to each character x^r Let us de- 
note Ti-^rj the set of vectors with other components equal to zero. Now, the 
equalities ( pp| ) are to be proved to end this section. 

Consider the decomposition of vector X(u, v) into a sum 

X(u,v) =^2X £in (u ) v), 
6,v 

where X^ G Wg,^. Then 

X^ v (g{u,v)) = xt, v (g)X^ v (u,v). 

Consider the sum 

E X^Gr'MM + B-uvC- vD)X(u, «)} = (21) 
see 

(which is equal to zero because of (|l^) ) . The action of g upon the braces in 
means that each u and v in the braces is transformed according to ([16]), ( 
i.e. u changes into Xi,-i{g) u i an( i v changes into X~i,-i{g) v - The equality 
gives thus 

uAX^i n +i(u, v) + BX^ n (u, v) - uvCX^ v+ 2{u, v) — vX^ + i v +i(u, v)D = 0. 

(22) 

Let us set u = in (^2|). Then v can take n different values Vj satisfying 
relation P(0, Vj) = 0. To these values Vj correspond d linearly independent 
vectors X^ v (0,Vj), and also d vectors X£+i lJ7 +i(0, Vj). Thus, the equalities 

BX^(0,Vj) = VjDX^ + i !V+ i(0,Vj) 
9 




that result from (22) give 



Analogously, one can as well obtain the rest of equalities ( poj ) . 

4 Connection to dimer model 

As has been demonstrated, the integrals of motion of the dynamical system of 
Section |l| and its reductions (if the even degrees of the transformation (0) are 
considered) are the coefficients a,j% of the vacuum curve @. These coefficients 
are determined up to a common factor, so they may be divided by aoo- As one 
can see, the resulting coefficients are those of the polynomial 



In other words, the determinant ( p3[ ) is an integral of motion for any u, v. 

Let us turn now to the model from section |3], that is to the model in 2+1- 
dimensional discrete space-time with periodic boundary conditions, and let the 
dimension d of the linear space corresponding to each vertex be equal to 1. 
Each of the "small" matrices A, B, C, D corresponding to the links will then 
be a single (depending on the link) number a,b,c or d. It is well known that 
the determinant of any N x N matrix is a sum of its matrix elements products 
corresponding in a certain way to the permutations of N objects, while each 
permutation decomposes into a product of the cyclic ones. In our situation, the 
cyclic permutations correspond to the non-sclfintersecting closed paths (con- 
tours) going along the arrows of the following diagram (Fig. |2|) (thus, general 
permutations correspond to the sets of non- intersecting paths). To each closed 
path corresponds the product of the weights ua, —uvc, —vd, b~ l on its links, 
and, to get right signs for the terms of which the determinant (^3| ) is made up, 
one should add a minus sign to each such product containing an even number 
of the factors b . 

Remark 1 Another way to obtain right signs is: to multiply each b by — 1 and 
then multiply each product corresponding to a closed path ( and containing any 
number ofb's) by — 1. 

It turns out that the determinant (j2jj) is connected with the statistical sum 
of the well known dimer model ||. Let us define the correspondence between 
the sets of paths and the dimer configurations as follows. Let the empty set 
of paths correspond to the "standard" dimer configuration, the dimers being 
placed on the "i?-links" (Fig. |J). For a non-empty set of paths, let us change 
the standard configuration along all the paths, replacing each dimer by a free 
link and vice versa. One can verify that this is a bijective correspondence. 




(23) 
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u a 



-uvc 
—vd 



ua 

-uvc 
—vd 



-0- 



—vd 



aa 



T 



-uvc 
—vd 



-uvc 



Figure 2: The ways along these arrows are connected with both the vacuum 
curve and the dimer model 



Figure 3: The standard dimer configuration 



The statistical sum being considered, let the weights —b (not corre- 
spond to the "£>-links", while to the other links correspond the unchanged 
weights ua, —vd, —uvc. Then one can see that the statistical sum, if multiplied 
by Ilover all links us cau the result the normalized statistical sum), 
consists of the same terms as the determinant ( p3| ) , up to different signs of some 
of them. Let us emphasize that the dimer model is, of course, inhomogeneous: 
the weights a, b, c, d are different for different links. 

Let us study these signs in detail. Note that the conditions of non-inter- 
secting and non-selfintersecting impose strong restrictions on the possible path 
configurations. Every closed path on the torus is homologically equivalent to a 
linear combination with integer coefficients of two basis cycles a and b whose 
intersection number is 1 (I use the boldface font for cycles, because the letters 
a,b . . . are already in use). If the torus is cut along a closed non-selfintersecting 
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path c not equivalent to zero, the result will be homcomorphic to the lateral 
surface of a cylinder (this follows, e.g., from S, chapter 1, section 3). Then the 
contour d going along a generatrix of the cylinder in a properly chosen direction 
has the intersection number 1 with the contour c. The intersection number 
being bilinear and integer-valued, we find that if the contour c is homologically 
equivalent to a sum la + mb, then I and m cannot have common divisors (not 
equal to ±1). Thus, the following lemma is valid. 

Lemma 1 Every closed non-selfintersecting path on the torus is homologically 
equivalent to a linear combination of the basis paths a and b with relatively 
prime integer coefficients. 

Now let us pass to the case of several contours on the torus. If they do not 
intersect, their intersection numbers equal (of course) and thus their homolog- 
ical classes must be proportional to one another. This together with Lemma [l] 
leads to the following lemma. 

Lemma 2 Several closed non-intersecting and non-selfintersecting paths going 
along the arrows on the torus, as in Fig. 0, are necessarily all homologically 
equivalent to one another. 

If two paths are homologically equivalent, then the terms of the same degrees 
in u and v correspond to them (one can see in Fig. ^| that the different ways 
round an "elementary square" yield the same degrees of u and v). Let the 
basis paths a and b yield the terms proportional to x = u ai v^ x ,y = u a2 v^ 2 
correspondingly (with the factors of proportionality not depending on u, v). 
According to Lemma [y the determinant ( p3|) and the statistical sum of the 
dimer model are polynomials in x, y. The following lemma sums up this section. 

Lemma 3 Let f(x,y) and s(x,y) be the determinant and the normalized 
statistical sum of the dimer model considered as functions of x and y. Then 

s(x, y) = \(-f{x, y) + f{-x, y) + f(x, -y) + f(-x, -y)) , (24) 

f{x,y) = ~(-s(x,y) + s(-x,y) + s(x,-y) + s(-x, -y)). (25) 
Proof. If the normalized statistical sum consists of the terms 
c jk x j y k = c jk (u a 'v^y(u a2 v^)\ 
then the determinant consists of the same terms multiplied by 

^ -^number of contours ^ ^■jg.c.d.(j.fe) /_ -^yk+j+k 

(here Remark [l] and Lemmas [l] and ^ are used) . This means that the signs of 
all the terms must be changed except where both numbers j and k are even. 
This is exactly what the formulae (p4|, E5I) do. The lemma is proved. 
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5 The discrete analog of Lax pair and a gener- 
alization of the dynamical system 

Now let us return from the reduction of Section p] to general matrices L = 
( A B" 
\C D / 
viewpoint. Denote 



Let us consider the evolution described in Section [l] from another 



A B 
C D 



This means that 



It follows from the equality ( pq ) that 



( A 








\ B 


D J 


(c 


5) 



(26) 



AA + CC = BB + DD, (27) 
AB + CD = 0, 
BA + DC = 0. 

These three equations are equivalent to the fact that the following equality holds 
for any complex u: 

-(A- uB)-\C - uD) = (uA + B)(uC + D)- 1 . (28) 

Vice versa, from (E8h follows 



L L = 



F 
F 



F being equal to both sides of (f27j), i.e. 



Mo I )<-')'■ 

It is clear that with any choice of F the matrix L belongs to the same equivalence 
class. The formula defines the same evolution in the space of these classes 
as it was in Section [jj with the agreement that the operators without a tilde 
correspond to the moment of time r, while those with a tilde correspond to the 
moment r + 1. 

The formula ( p8| ) suggests the following generalization. Let, from now on, 
A(u) and B(u) be matrices depending polynomially on u: 

A( U )=A +A + --- + ^m^ roA , (29) 
B(u)=B +B 1 + ... + B mB u mB . (30) 
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We will look for matrices A(u), B{u) — the matrix polynomials of the same de- 
grees rriA and ms in u — that satisfy, for any u, the equation 

B(u)- 1 A(u)=A(u)B(u)- 1 . (31) 

The relation (^l|) provides what is called a discrete analog of the Lax L, A-pair, 
which means here that the operators A{u)B{u)~ l and A(m)_B(u) _1 (which are 
playing the role of L of the pair) are "isospectral deformations" of one another: 

A(u)B{u)- 1 = A^A^Biu)- 1 !^)- 1 . 



Let v be an eigenvalue of both sides of (|3l|). Let Y(u, v) be the corresponding 
eigenvector normalized, as in Section ||, so that its last coordinate equals unity, 
and let X(u,v) be the vector proportional to J3(u) -1 Y(u, v) and normalized 
in the same way. One can verify that this may be described by the following 
formula (h(u,v) being a scalar factor): 

( \ X(u,v) = h(u,v) ( I \ ®Y(u,v), (32) 

which is in obvious analogy to (^). The divisor equivalence is 

mD u + D x ~ D v + Dy, (33) 

D u and D v being pole divisors of the functions u and v, m = max^i^B). 
For a given u, the eigenvalues v come from the equation 

P(u, v) = det(A{u) - vB(u)) = 0. 

It defines an algebraic curve T — "generalized vacuum curve". Let us calculate 
the genus g of the curve T. First, we need to know the number of branch points 
of the projection 

(u,v) — >u (34) 

of the curve T onto the complex plane. 
Consider P(u,v) as a polynomial in v. 

P(u, v) — a Q (u) + a 1 (u) + . . . + a n (u)v n . (35) 

One can verify that aj (u) has a degree 

degcij(u) = (n- j)m A +j B - (36) 

From this one can deduce that the discriminant of P(u, v) considered as a poly- 
nomial in v is a polynomial of degree 

b = (rriA + rn,B)n{n — 1) 
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in u. The mapping ( p4[ ) being n-sheeted and the number of branch points 
equalling b, one obtains from the Riemann — Hurwitz formula that 

/ Y s f m A +m B \ , , 

g=(n-l)l n-lj. (37) 

So, the following construction has been described. Given two polynomial 
matrix functions A(u) and -B(u), one considers the meromorphic matrix func- 
tion A(w)_B(u)~ 1 (or else B(u)~ 1 A(u)), and from this function the algebro- 
geometrical objects arise: the generalized vacuum curve T and the linear equiv- 
alence class of the pole divisor Dy (or, respectively, Dx) of the eigenvectors of 
the mentioned meromorphic matrix function. Instead of the pair (A(u), B(u)), 
it is sufficient to indicate its equivalence class with respect to gauge transfor- 
mations 

A(u) -> GA(u)H, B(u) -> GB(u)H; (38) 

instead of the function A{u)B(u)~ 1 , its equivalence class with respect to trans- 
formations 

A{u)B{u)- 1 -> GA{u)B(u)~ 1 G~ 1 

will suffice. Then it turns out that the correspondence between such equivalence 
classes (either of the pairs (A(u), B(u)) or the functions A(u)i?(-u)~ 1 ) and the 
abovementioned algebro-geometrical objects is a birational isomorphism, the 
divisors Dx and Dy being of degree g + n—1, as in Section |[ 

The easiest way to show this is to start from a given curve T defined by the 
equation 

n (n-j)m A +jm B 

P(u, v) = ^2 ^2 a jk v j u k = 

(compare with (|35|, ^)) and a divisor Dx in it of degree g + n — 1. The number 
of coefficients djk minus one common factor equals 

(n + l)(^±^n + l)-l. (39) 

The linear equivalence class of divisor Dx is defined, as is known, by g param- 



eters. Adding up the expressions (B9j) and (37), one gets the total of 



(m A + m B )n 2 + l (40) 

parameters. 

Then, the gauge equivalence class of the pair (A(u), B{u)) is constructed out 
of relation (]3^). To give more details, one must at first choose a divisor Dy 
satisfying the equivalence ([33|). Then the poles and zeros of the function h(u, v) 
are determined. For X(u,v) and Y(u,v) one must take columns consisting 
each of n linearly independent meromorphic functions with corresponding pole 
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divisors. The arbitrariness in these constructions leads exactly to the fact that 



A(u) and B(u) are determined up to a transformation (38). 

The pair (A(u), B(u)), up to a scalar common factor, is determined by + 
hib + 2)?2 2 — 1 parameters (see ( p9| , p0|)). In taking the gauge equivalence class, 
the number of parameters is reduced by 2(n 2 — 1). The result is again (fIo|). 
This means that, indeed, to a generic pair (A(u), B(u)) corresponds a divisor 
Dx of degree g + n — 1 and the correspondence 

(gauge equivalence class of the pair (A(u),B(u))) < — ► (T, class of Dx) 

is a birational isomorphism. 

Now let us recall that Y(u, v) was defined as an eigenvector of the operator 
A(u)_B(u) -1 , while X(u,v), as is easily seen, is an eigenvector of B(u)~ 1 A(u). 



The relation (31) means that for the pair (A(u), B(u)) its vector X(u,v) is 



nothing else than Y(u,v), i.e. the equivalence holds 

D x ~D x + (mD u -D v ). (41) 

Now, assuming that if a quantity without a tilde corresponds to the moment 
of time t then that with a tilde corresponds to t + 1, one comes to a conclusion 
that to the adding of unity to the time corresponds a constant shift ( f4l| ) in the 
Jacobian of the curve Y. Thus, the dynamics of the system in this section, as 
well as in Section 0, linearizes. 



6 Discussion 

In this paper I study a dynamical system in discrete time, i.e. a mapping and its 
iterations, acting on finite sets ofuxn matrices. The system appears in several 
modifications, on which depends the number of matrices as well as the additional 
conditions that may be imposed on them. The "law of motion" is formulated 
in a rather simple way, and a large number of "integrals of motion" turn out 
to exist and be the coefficients of the "vacuum curve" — the object coming from 
the theory of the quantum Yang — Baxter equation. If the motion in the system 
is considered up to a "gauge transformation" , the system is integrable in the 
sense that there exists a birational isomorphism between the "phase space" and 
the set of pairs (a vacuum curve, an element of its Picard group), so that in 
the process of "motion" the vacuum curve doesn't change, while the element 
of its Picard group depends on the time linearly. Thus, the Cauchy problem is 
solved through the following scheme: the initial point in the phase space — ► 
the vacuum curve and the element of its Picard group at the initial moment of 
time — > the same at the moment t — ► the element of the phase space at the 
moment r. 

Connections with statistical physics are exposed in Section where a special 
reduction of the model is considered. Note that there exists one more dynamical 
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system in 2+1-dimensional discrete space-time that is also connected with sta- 
tistical physics and seems to be completely integrable. This system is as follows. 
Consider the inhomogeneous Ising model on the triangle lattice. Imagine this 
lattice as consisting of triangles of the form A and perform for each of them 
the "triangle — star" transformation 0: 

Thus, the hexagonal lattice appears, which now may be imagined as made up of 
its parts of the form "V . So, let us perform the "star — triangle" transformation 

for each of them. One step of the evolution is over. It would be of interest to 
reveal possible connections of this model with the model on the square lattice 
from Section ||. 

Another interesting problem still unsolved: to describe the evolution of the 
system in Section [I] in full, not up to gauge equivalence. 
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